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NATIONAL AD7IS0RY COMMITTEE FOE AERONAUTICS 
TECHNICAL NOTE NO. 1171 

ON SUBSONIC COMPRESSIBLE FLOWS BY A METHOD OF CORRESPONDENCE 
II - APPLICATION OF METHODS TO STUDIES OF FLOW WITH CIRCULATION 

ABOUT A CIRCULAR CYLINDER 
By Shepard Bartnoff and Abe Gelbart 

SUMMARY 


A general method for studying the flow of a compressible fluid 
around a closed body has been discussed previously in part I of this 
report. In the present paper application iB made to the specific case 
in which the linearized equation of state is used. For a given incom- 
pressible flow around a specific profile, a corresponding compressible 
flow is found. The flow at infinity remains unchanged. Detailed stud- 
ies are made of the flow with circulation around a unit circle, and 
velocity distributions are found for a wide range of Mach number and 
angle of attack. Comparisons are made with other methods. 


INTRODUCTION 


■ The present report 1 b the continuation of a previous report by 
Gelbart (reference l) in which a general method for studying the flow 
of a compressible fluid around a closed body is discussed. The method 
is based on finding compressible flows that correspond to given incom- 
pressible flows. 

Since the compressible complex potential in the general case is 
not an analytic function, the ordinary theory of analytic functions of 
a complex variable is not applicable. However, in the hodograph plane 
(where the variables are the velocity magnitude and the direction of the 
flow) the complex potential is a function of the type studied by Bers •. 
and Gelbart (references 2 and 3) and termed by them "sigma-monogenic . " 

In the present report, the condition under which the differential 
equations of a compressible flow in the hodograph plane become Cauchy— 
Riemann equations is used. This occurs when the linearized equation of 
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state is vised (7 = —l) . The complex potential in the holograph variahleB 
is then an analytic function and the theory of functions of a complex 
•variable is applicable. The linearized equation of state is used 
throughout in this report . 

Emphasis is placed on thA compressible flow with circulation around 
a unit circle. The correspondence function is chosen so as to yield 
this flow to a very close approximation even for fairly high free- -stream 
Mach numbers. The distortion of the body is very slight, but, as is to 
be expected, increases for very high free-stream Mach numbers . The re- 
sulting body shapes are studied for the entire range of possible free- 
stream Mach number, which for subsonic flows is from 0 to 1. Resulting 
body shapes are also studied for different angles of attack. Finally, 
velocity distributions are computed around the slightly distorted unit 
circle for different free— stream Mach numbers and for different angles 
of attack. 

Some comparisons are made with the results of Tsien and Eers. 

Bers’ and Tsien' s formulas turn out to be special cases of the main for- 
mula derived in this report. One advantage of the present method is 
that it yields flows with circulation about nonsymmetric closed profiles. 

This work was sponsored by and conducted with the financial assist- 
ance of the National Advisory Committee for Aeronautics . This report 
was submitted in July 1945. 


GENERAL FLOW THEORY 


Several of the more commonly used equations governing the behavior of 
fluids will be mentioned here without proof. These are relations involv- 
ing the quantities: velocity q, pressure p, density o, ratio of the 

specific heat at constant pressure to the specific heat at constant volume 
7 , and the velocity of~sound a. Subscripts of zero (i.e., p Q , p Q , a Q ) 

refer- to values of the respective quantities at a stagnation point (q - 0 ) . 

The first fundamental relation is Bernoulli's equation, which may be 
written in differential form 


4 



0 


( 1 ) 


Another fundamental relation is the isentropic relation (adiabatic 
equation of state) 
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P = Po 

The velocity of sound is given by 


f-2-X 


\fioJ 


a 2 = 4P _ TV 


Hence 


dp 


a o 2 = 7 


Po 

Po 


The units of density and velocity will he so chosen that 


P 0 = 1 


a o " 1 


The quantities p and q will bhen he dimensionless , 
The following relations are easily derived. 


( 2 ) 


( 3 ) 


W 


( 5 ) 


! = i - 1 (7 - i) a 2 

- (6) 

- i 

1 - | (7 - 1) <1 £ 7-1 

(T) 

r . 7 

,[1-| (7-1) ^J 7 " 1 

(8) 


The Mach number M 1 b defined by 


M ** — 
a 


( 9 ) 


The symbols and q^, refer to the values of these respective 

quantities in the free stream. 
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The velocity potential 0 end. the stream function vj/ for an incom- 
pressible flew satisfy the differential equations. 



the subscripts x and y denoting partial differentiation with respect 
to these variables. These equations are the Cauchy— Hi emaan equations, so 
that 


fj = '0 + iij/ 


( 11 ) 


is an analytic function of the complex variable z » x + iy. « 

If instead of using the independent variables x and y in the 
physical plane, the independent variables q and 9 in the hodograph 
plane are used, where q Is the magnitude of the velocity and 9 is 
its direction angle, the differential equations became 


0 9 " <1 ^ 

- - !+• 


; 


( 12 ) 


By making the change of variable 


q 



the system (12) becomes 


09 = ^q 
02 “ - ^e 


(137 


( 14 ) 


which are again the Cauchy— jRiemann equations . Thus, in the hodograph 
plane the incompressible complex potent ialL-is an analytic function of 
the complex variable w = 6 + iq. 
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The physical and hodograph variables for an incompressible flow are 
connected by the relation, 


dz = e lw dQ 


(15) 


The analogous differential equations for a compressible flow in the 
physical plane are 





1 

p ^ 


and in the hodograph plane are 


(16) 




Sk- 

p 





1 - M 2 

pq. 



(17) 


System (16) is nonlinear, and no systematic mathematical treatment for 
such equations exists. However, methods for a systematic treatment of 
the solutions of equations (17) have been developed in reference 1. 

The equation analogous to equation (15) for a compressible flow is 


dz 


■ie 
S — 

q. 


(<¥ + i d\l/^) 


(IS) 


Equations (17) may be symmetrized by the change of variable 



( 19 ) 


Note that q = 0 . 

05 

Equations (17) then become 
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.yr~_Mg 




yi - m s 


^e 


( 20 ) 


Xf = constant, equations (20) reduce to the Cauchy-Eiemann 


equations . When 


Jx.-.ut 


constant, the constant can he evaluated hy 


choosing q * 0 . Here M 0 = 


M(q) 


_ = 0 ariri 

q=o 


JT- m £ 


q=o 


Po 


“ 1 


(21) 


This implies that y = ~ 1. For when M is replaced hy q/a, it 
follows from equation (21) that 


a 2 — q 2 — a 2 p 2 * 0 


Then, hy differentiation. 


2a g (1 - p 2 ) - 2q - 2a 2 p ^ = ° 


From relation (6), 


2a ~ - -(7 - lh 
dq 


and from relations (?) and (6), 

d£ 

dq 


— — a (i - i e ) 1 




7_nJL n £ 


ill 




q.p 


( 22 ) 


(23) 



a 
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Substitution of these values into equation (23) gives 


[-(7 - 1) - 2][(1 - p a )q = 0 

[- 1(7 - 1 ) - 2 ] -■ 0 

77 = -1 (24) 


With this value of 7, the equation of state (2) becomes 


-1 

P = Po P 


(25) 


and since the density is proportional to the inverse of the volume, the 
volume becomes a linear function of the pressure. This linearized equa- 
tion of state will be used throughout the remainder of this paper. 

It should be emphasized that no actual gas satisfies the pressure- 
density relationship described by 'equation (25), For actual gases, 7 
lies between 1 and 1.6, having a value of 1.4 for air. The simplifica- 
tion resulting when 7 = —1, and the fact that the theory of analytic 
functions of a complex variable then becomes applicable, led Chaplygin 
and later von Barman, Tsien, and others to Btudy the properties of such 
a fictitious gas (references 4 to 11) . 

Von Barman (reference 7)j however, showed that using the 
value 7 “ — 1 is equivalent to replacing the curve of pressure against 
reciprocal density in the adiabatic case by a tangent line to that curve. 
At and near the point of tangency, the fictitious gas approximates the 
behavior of the actual gas. The study of the behavior of the fictitious 
gas is Justified by the insight such a study gives into the solution of 
the actual problem for subsonic flows . 

When 7 * — 1, equations (6) tc (8) become 


a 2 =* 1 + q 2 
1 

p = — 

yiTF 

p * p o y 1 + q 2 


(26) 

(27) 


(28) 
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and it follows directly that 

M 

q. = •- ===■— 

J\ - M 2 

Equations (20) reduce to the Cauchy- -Riemann equations, 

(*q - ~ * e 

Thus, Qa <f + i\j/ is an analytic function of v » 9 + iq. 
From equations. (19), (21), and (27) it follows that 



/ \ 
■ 1o « r — 7==a=) 

\ Wl + q a/ 


where 

K s t ^/l *t~ Sco 
1 

03 

Hence, 

e ? - ■& 

1 +Av/l + q 2 

From this equation, 

1 e K _ 
q 2e^ 2K 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


Also, from equations (33) and (27), 
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a + - 1 


Pa 


- X 


e^ q. 


and by use of equation (3^)> 


1 K t e u 
pq. 2e^ 2K 


(35) 


Upon substituting the values of ^ and — from equations (3^) and 

» * * 

(35) into equation (l8) # it follows that 


dz 


ie 

C 

2 

»ie 


;(J - x)^ +i (5 + f) 4 *] 


-4 (d0 + id\|/) - ~ (d0 - id^) 

04: K 


| e 1 (e + ia) (d$ + id\J/> - i e i (e ~ 1$) (d0 - id\|/) 


E ^iw . 1 — iw 

- e dQ- — e dQ 


and 



(36) 


Since Q = $ + ity is an analytic function of the complex variable 
w = 0 + i§", the mapping 


w =.-i 


log 


a 

K rj®(£) 


( 37 ) 


defines a complex potential G(£) of a compressible flow in the £ — 
plane,, where f(£) and G(£) are analytic functions of the con® lex 
variable £. The region of regularity of f(£) and G(£) will be 
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dealt with later. The function. f(£) will "be called the correspondence 
function. 

From equation (37), with primes denoting differentiation with re- 
spect to £ p 


and 


iw _ 2 f * LU 

“ K G*(£) 

(38) 

-Iw KG'(f) 
“ 2 f«(£) 

(39) 


After substitution of equations (3$) and (39) into equation ( 36 ), it 
follows that 


so that 




5 / 4£ 


-i / l f, «> 


= f«> -£ j' 


(40) 


If G(£) represents the flow around a given closed body in the £— 
plane, equation (40) will map that flow into the z-plane. The corre- 
spondence function f (£) will be chosen in such a way that the flow in 
the z-plane will be essentially around the same profile as in the £ — 
plane . 

The compressible-flow velocity in the z— plane may be determined 
from equation (39) • For, 


iw _ —i (e+iq) _ K G« ( t) 
2 f 1 (£) 


-ie q K 
e e = — 
2 


0L£SI 

fHO 


(41) 


By taking absolute values, 
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£ K 


K 

Q *(0 

2 

«(£) 


and making use of equation (33)* 




1 +yi + <r 


1 

2 




it follows that 


k 

& T (D 



JJILL 


4- 

5!itl 

f*(S) 

2 


(42) 


Stagnation points occur when q = 0. Thus, the stagnation points 
are located where 


<**(£) - 0 


(43) 


It is desirable that the correspondence should not alter the flow at 

Infinity. At infinity, 6=0 and = 1. An examination of equation 
(4l) shows that the condition for the flow to remain unchanged at Infinity 
is 


lim KG f (0 _ 

£-»•» sfTJTT- 1 


(44) 


Since G* (£ ) Is regular and different from zero at infinity and if 
g '-f is " fco ©xist and be different from zero, the most general 

form that f 1 (£ ) can have is 


f*(£) 


V 


bn 


■ n* 


t G 4 0 


(45) 


n=o 


where 


l = re 


IX 


Then 


(46) 
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w f 'b 

f (£) ■ fc-1 + lo £ + t>! log £ - > ,% (47) 

fe £ 

Also 

lim f»(£),i> 0 

t — >0* 

and from equation (44) 

"o" t ^„| «'••(£) (48) 


Thus, the condition of the flow at infinity in general predetermines the 
value of b Q . 


It is desirable to obtain the integrand in the right-hand side of 


equation (40) as a power series in 
Set 




n=o 



n=o 



(49) 


Then 



Equating coefficients of like powers of 



n 
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1 - *o B 0 


0 — bj_ Bo + bo Bx 
0 — b 2 B 0 + bi Bj. + bo B 2 


0 0 B n B o + B i + • • • + Bj^ + b 0 B n 


These equations yield the iteration formulas 
b ° = £ 

Bn = — ^(Bn B o + Bn-x B i + . . . + b 2 Ti^z + bx Bj^) 
D o 

where n > O, or 


n n+x 

*0 


Bo 

0 

0 

0 

• • • 0 

Q 

0 

0 

1 

Bx 

Bo 

0 

0 

9 9 9 0 

0 

0 

0 

0 

b 2 

Bi 

Bq 

0 

... 0 

0 

0 

0 

0 

* 3 
• • 

b 2 

• • 

Bx 
• • 

Bo 

• • 

... 0 
• • • • 

0 

9 

0 

• 

0 

• • 

0 

• 

• • 
Bn— 3 

• * 
Bn— i 

• • 

Bn-e 

9 9 

Bn-e 

• • • • 
* . . Bj_ 

• 

Bo 

• 

0 

• « 
0 

• 

0 

Bn-a 

Bn-3 

Bji — 4 

Bn— e 

« • • l3g 

Bx 

B 0 

0 

0 

^n-i 

^n-a 

^n — 3 

^n—t 

... hj 

b 2 

Bx 


0 

Bn 

Bn— x 

Bn- a 

Bn — 3 

• • « b 4 

b 3 

B a 

Bi 

0 


Then set 
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_i / ISLU I I 
f*(0 


00 

n«=o 

CO 

» Cl I + Co log J + ^ C n Aj 

n=i ^ 


(53) 


Equation (40) then 100011106 


= b_! + i 0 £ + b! log£ - ( f is) 

'£=1 t ' 

OO 

+ C— 1 ^ + Co log £ + Cn-~Q 

n=»I * 


(54) 


i X. — ^ \ 

When £ is replaced, hy re and £ by re , equation (54) beoamos 

z =■ b— ! + b 0 re 1 *" + b x log r + bi(iX) 

OO 


/ V ^n+i 1 -inX\+ C — 1 

" vZ, “5” 7* 6 ) 


re ^ + C c log r 


n=i 


c 0 (— ix) + y 1 C n ^ 


F e 


inX 


(55) 


n=i 


If a dosed contour in the £-plane is to be mapped by equation (55) 
into a closed contour in the z-plane, 


bi - C 0 = 0 


(56) 


This predetermines the value of bi. 
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Set 


N =* t—! + bj. log r + C 0 log r 
Under the condition of equation ( 56 ) 

( f 

n=a 


+ 


Oo 

r. -il 
C_i re + 

n=x 




ini 


(57) 


(58) 


TEE FLOW AROUND A CIRCLE 


The example of a compressible flow that will he treated here in 
detail is the flow with circulation around a circle. This example haB 
been chosen because of its fundamental importance. 

In determining the correspondence that will give the compressible 
flow around a unit circle, the complex potential function chosen is that 
of an incompressible flow with circulation around a circle of radius R 
in the £— plane, namely. 


OU) = *,(£ -glog| 


(59) 


where F is the circulation, 
(equation (43)). 


G ‘(0 


Stagnation points occur where 

m 1 


(**(£) -.0 
(60) 


il 


If 


£ = Re 


> 
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G* = q _ IE 1 e ^ . 5 e 
2n E 


-ieA 


" 4 »~ if 5 003 X - h | 3ln X - 4„<»<r 8X + sin 2X 


= Q. (l **• cos 2A. ) — — 


a, e aln x - 1 (s ^ 2X - £ | «« x 


•) 


“ sln X ( 2 V aln x ~ |) + 1 °°o x sin X “ ^ 5) 


= 0 


Hence 


IX 


( s % 


sin X H 

00 2k E/ 


I G T | = 2q sin X - i 
« 2k B 


Let there he stagnation points at X = a, X ® l80° - a. Then, 


2q sin a E- — => 0 

00 2k R 


— = 2q. E sin a 

dJl OO 


(61) 


Thus., fixing the circulation fixes the stagnation points and hence 
the angle of attack. 

Equation (48) is used to fix the velocity at infinity. Since 


lim G * (t) * q 

— > 00 c 

_ *1® 


(62) 


b ° " 2 


t(yr7 


CL 2 ^ 1 ) 

CD / 


(63) 
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The equation of a unit circle is £ = e If equation (58] is to 
approximate such a circle, the sum of the coefficients of the e —terms 
should equal 1 and the coefficients of the other terms should approximate 
zero. The first step in achieving this end is to set the coefficients of 
all the e _in?k '-terms, n > 2, equal to zero; that is. 


bn - 0, n > 2 

* 

The correspondence function f(£) then reduces to 


(6k) 


and 


f(£) = b~i + £ + hi log £ - b ; 


f *(£) = b 0 + b x i + h 2 \ 

t t 


(65) 

( 66 ) 


Also, equation (58) beoomeB 


w , . -j iX. b 2 . —ix, 
z = N + b Q Re ~ — e 


+ C— 1 Re ^ + Ci g e ^ Cn 


1 ini 

i * 6 


(67) 


n=2 


il 


The sum of the coefficients of the e — terms may be made equal to 1 
by fixing R so that 


b 0 R + Ci J = 1 


R 


( 68 ) 


Also the e term may be eliminated entirely by fixing b 2 so that 

b 2 


R 


+ C— >1 R =* 0 


(69) 


Then equation (67) becomes 


„ ^ il ^ V n 1 „ inX 
z — N + e + y Cn — g e 


R 


( 70 ) 


n=2 
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In this equation N is a translation constant that does not affect 
the shape of the tody. The Cn coefficients are shown to be small and 
therefore the resulting body shape is essentially a circle . 

It will be shown later that all the odd C a coefficients are real 
and the even C& coefficients are pure imaginary. Therefore 


z = x + iy 

c c 

«* cos X sin 2X + C cos 3X — sin 4X + C_ cos 51 + • • • 

i 3 i 6 

+ i^sin X + cos 2X +• Cjj sin 31 + ^ cos 4X +- sin 51 

+ . . . \ + N Cn) 

where the constant N is chosen tc make the absolute value of- z when 
X = 90° equal to the absolute value of z when X = 270°, 

N =• C 2 — C 4 + C 6 - C e + . . . (72) 


and N is pure imaginary. The resulting body is symmetric with respect 
to the y-axis. 


After the circulation is fixed to give a desired angle of-attack by 
equation (6l), expression (59) for the complex potential . G(£) reduces to 


and 


Then 


G(t) = (i - i 2R sin a log | 

<**(£) = q w (l - i 2 sin a - | 5!) 

2 

[ a* (z)] - A x - iA £ — - A 3 - — + iA 2 — + A x 


E 4 

r 


(73) 

( 74 ) 

(75) 


whore 
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Ai = q. 2 

j* ^ CO 

A 2 = 4q. 2 sin a 

CO 

A = q. 2 (4 sin 2 a + 2) 

3 co 

and 

C— i ** — ^ A x Bo 

C 0 = - | (A x Ii - iA^oE) 

Ci = ^ (Aila - iA^iB - AjBqH 2 ) 

c = ■— 1 — (a x B- - iAaBgR - A,B X B 2 + 1A^ 0 J? ) 

4.2R 2 3 

C = .i.- -(AA r* iA-JB-B - A-BaB 2 + iA^ x B 3 + A X B 0 B 4 ) 
3 4.3B 3 

c 4 ~ ^ CA X B S - JLA2B4 B — Aj^RjJP + iA 2 B 2 B + A x BiB ) 


( 76 ) 


(77) 


C n = (AiB^j. - iAsB n B - A L / + iAsB^^R 3 + AiB^R 4 ) 

4.rdC 


Also, since b n = 0, n > 2, the relations for B n as given in 
equations (51 ) to (53) become 



Bx- - - h- (b x B 0 ) 
D o 


B2.* — tb2B 0 + b x Bi) 
B 3 = — ^ (b2B x + b x B 2 ) 


( 78 ) 


B 


n 



(baB jj _ 2 + biBQ-i), 


n> 1 
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In determinant form, 


B r 


n+i 


to 0 
*1 b 


0 

0 

• 

0 

0 

0 

0 


0 

0 


b2 b^ bo 


b 2 

0 


bx 

bg 


0 0 
0 0 
0 0 
0 0 


0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 


0 b 2 bj. bo- 0 

0 0 bg bj_ b c 

0 0 0 b 2 bi 


1 

0 

0 

0 

0 

0 

0 

0 


0 0 


bg b x 0 


(79) 


.Also Bn is given by the summation 


Bn 


-T 


r=o 


(-l) n+r (n - r)t bx 


n- 2 r 


r! (n — 2r) l 


bo 


n-r+i 


(80) 


Thus, the first few Bjj-texms are 


B ° = b 0 ' 


Bx = 


bi 

2 


Br 


b 3 
°o 




(81) 


The first two C n — terms as given by equations (76) and (77) axe 
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C— i = — T7 Ai 


C 0 = + i ^ 

^o 2 ^ 


Equation ( 69 ) is nov used to determine Tig: 


- = 2 - + C-iE *= 0 
E 


5a -i a x A 

E 4 X b 


*bo - — 


o 

A^B 2 

4b„ 


Equation ( 56 } is used to determine bj.: 


"bjL. Co 36 0 


*1 


_ A l^i- _ i ^2® = 0 

4b 0 2 4b 0 


4bQ^bi i — A^b-^ — iAgtgE 


The quantity hi must he pure imaginary . Therefore, 

hj. = -hi 

4b ( - ) 2 bi + A^bx * iAghoE 
A a b E 

hj. = i p 

kb 0 + A x 


( 82 ) 


(83) 


(84) 
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If in like manner equation (68) is used to evaluate R, the result- 
ing complicated expression may he reduced to 

R - ^ (85) 

1 — M s ein a a 

It is interesting to note that all odd B^— terms and all even Cq— 
terms are pure imaginary, whereas all even B a -terma and all odd 0 n — 
terms are real; also, that when a = 0°, corresponding to a zero angle of 

attack, R « b . by =« 0, and tr + -8 » 1. 

1 O' x * o JJS 

Equation (42) is used to determine the velocity^ on the body corre- 
sponding to any value of The velocity is limited by the restriction 

Sliii % k (86) 

*•<£> 

When equation (6l) is substituted into equation (60), 

2 

<*•(£> " 4»(1 “ 21 B Bin ai-^) 

• q^l - 21 sin a e” 1X - e~ 2lX ^ 

= q 1 — 2 sin a sin \ — cos 2A 

°°L 

+ i(— 2 sin a cos \ + sin 21 ) ( 87 ) 

Equation (66) can also be rewritten 

f(£) - »„ + ^ <T lX + || e" 2U 
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The maximum value of 




f’(£) 


occurs when >. = 90 ° • For this value 


of G-* (£) = 2q. (1 — sin a,) and f*(£) = "bo + ~~ ~ -§• Hence, by 

s 00 iE B • 


restriction (86), 


2q. m (1 - sin a) 


< 2 


(89) 


bo + 


or *>£ 


iE E 

By substituting., from equatf.ons (83) and (8h) 

A p b 


q. (1 — sin a) < b a + 


?vo 


Ax 


(90) 


4Tj 0 £ + Ax toe 


By -substituting from equations ( 63 ) and ( 76 ) 

<3. m C 1 ~ slu < | (v^t + l) 

b /r < +1 ) , 


(V 1 + 


2c sin a ; 

4f./l + <T m + Z ) + <3 - 2 ® 2 (^ // 1 + ^a> + 3 ) 


(91) 


1 — sin a< 


1 + 1 


CO 

Z 

i 05 


sin a 


a 1 + q/L 

*00 CO 


and using equation ( 29 ), 


M < 1 + M (1 + M m ) sin a 

CO 00 oo 


( 92 ) 


Hence, for a zero angle of attack, M^, may vary between 0 and 1. 
For a given however, equation ( 92 ) limits the permissible angle of 

attack. 
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COMPARISON WITH TSIEN METHOD 


The equation derived fey Taien, reference 5, in the notation of the 
present paper, is 

z - ££•(£>]* H (93) 

•J 

Thus it Is seen that Tsien's equation is a special case of that used in 
the present report where 


*<£) 


Tsien's equation alters the velocity of the flow at Infinity, or the free— 
stream velocity, hut this is relatively unimportant since the amount of 
the change Is readily determined. However, since Talon's equation has no 
hi term, it cannot he used to study flows with circulation. In Tsien's 
equation r is therefore zero and G is 

o<£> - «»(? ■’■f) W 

Then, 


and 


In order to have Tsien's equation give as nearly circular a body as possi- 
ble, let 




V 


i 


2 


Z - t - i q s 1 - SihSf + 

s “ 5 n lai 3 


= Re 1 * 1 - ~ q 2 Re - J q e Re 

4 co ° x ~ 


IX 


2 co 


q 8 Re 


is X 


12 <» 




1 


(95) 
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Then, 


il 1 2 T3 — i^» ,3. 2 n i3^> 

z = e - 7 r Ee + — - q~ Se 

4 £» 12 0= 


= ( x ~ t ^ R ) cos x + cos 3X 

+ i j^l + | sin \ + i E sin 3x"| 


(96) 


Syracuse University, 

Syracuse, New York, July 1, 


19 ^ 5 . 
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TABLE I.- CONSTAHTS OT THE TBABSFQBMATIOH FOBCTIOff FOB VAHIOCS MACH BOMBEBS [a - -eO<3 


i 

0.1 

mm 

0.3 

o .4 

BOB 

0.6 

e 

0.100504 

0.204124 

0 . 3 i 44 s 6 

0.436436 

_ _ _ . _ 

0.577351 

0.750000 

E3 

1.002521 

1.010309 

i. 024 i 44 

1.045546 

1.077352 

1.125000 

B 

-t <0.003450) 

- 1 < 0 .oi 4173 ) . 

-i(o. 033399 ) 

-iC0,0636l8) 

-1(0.109575) 

- 1 ( 0 . 180760 ) 

B 

-0.00253s 

- .010623 

.025864 

.051706 

.095270 

.172420 

s 

1.003695 

1.015059 

1.035041 

1.065433 

1.109S07 

1.174459 

B 

-1(0.0017) 

-1(0.0069) 

-1(0.0157) 

- 1 ( 0 . 0283 ) 

- 1 ( 0 . 0458 ) 

- 1 ( 0 . 0722 ) 

°3 

0.000s 

0.0037 

0.007s 

0.0137 ' 

0.0212 

0.0308 

°4 


1(0.0003) 

1(0.0002) 

1(0.0005) 

-1(0.0031) 

-1(0.0067) 

°5 




- 0.0001 

0.0007 

0.0016 

B 





-1(0.0002) 

1(0.0006) 

°7 






-0.0007 

H 

-1(0.0017) 

-1(0.0072) 

-1(0.0159) 

— 1(0.0288) 

- 1 ( 0 . 0429 ) 

-1(0.0655) 


TABLE H.- COHSTAHTB OF THE TSABSFCTOttTICM TDUCTIOH 
K® VABIOOS AJC-I2S OF ATTACK [M^- 0.3) 


Si 

O 

O 

- 10 ° 

- 20 ° 


0.314486 

0 . 314486 

0.314486 

m 

1.024144 

1 . 024 i 44 

i. 024 i 44 

m 

0.0000 



b 2 

-0.025322 

0.025460 

- 0.025864 

E 

i. 024 i 44 

1.026931 

1. 035041 

°2 

0.0000 

1 ( 0 . 0080 ) 

-1(0.0157) 

°3 

0.0069 

0,0078 



0.0078 

°4 

0.0000 

s 1 

i (0.0002) 

°5 

0.0001 

Bl 


H 

0.0000 

1 ( 0 . 0078 ) 

•1(0.0159) 







































































































TABIZ m.- COBSTAJTOB 07 THE 
yarxoos AKnas or , 





0° 

-5° 

-10° 

0.980196 

0.9S0196 

0.980196 

1. 20011*0 

1.200l4o 

i.apoi4o 

0.0000 

1 

-1(0.072041) 

-1(0.145138) 

-0.285269 

-0.290427 

-0.296900 

1.200140 

1.204624 

1.218138 

0.0000 ■ 

-1(0.0249) 

-1(0,0496) 






0.055b 


0.0000 


n rw£ 


0,0000 


0.0007 


°8 

0,0000 

C 9 

0.0001 

T 

0.0000 


0.054s 


-i(o.oo 4 i) 


n rwvrz 


-1(0.0007) 


0.0006 


- 1 ( 0 . 0001 ) 


0.0523 


-1(0.0080) 


a nnh£ 

V • Ul/TV 


-1(0.0013) 


0.0004 


- 1 ( 0 . 0002 ) 


A 

-90 0 


V 


* 

*J ~~ 


l B J 

0,0000 I -0.9992 I 0.9992 I -90.0 0 


“75° 

0.2574 

-0.9656 

0.9993 

-75-1° 

-6o° 

0.4977 

-0.8668 

0.9995 

1 

-60.1° 

-45° 

0.7048 

-0.7080 

0.9990 

-45.1° 

-30 0 

O.S645 

-O.50I7 

0.9995 

-30.1° 

0 

S' 

1 

O.9656 

-0.2610 

1.0002 

-15. l q 

0° 

VQQQg 

_n nmli 
" * 

1 AAA® 

* • WWW 

_ n oO 

— V«fr 

15° 

0.9674 

0.2562 

1,0007 

i4.g° 

30° 

' 

0.8675 

0.4982 

i.ooo4 

29.9° 

45° 

0.7082 

0.7060 

1.0000 

45.0° 

6o° 

0.5007 

0.8652 

0.9996 

59.9° 

75° 

0.2591 

0.9651 

0.9993 

75.0° 

90° 

0.0000 

0.9992 

0.9992 

90.0° 






































































































t 


V 


3«bte T-— irragY5®0K OF pgryme ygoii GIHOEB 
0.2j tt *> -20°] 


-0.9963 


TABU 71.- YEOCirr KIB9B35D*rOB OB IBOTIII [h^- 0 . 3 } a - 0°] 



-75° 0.2524 -0.9643 


-45«> O.6976 -0.7172 


-30° 0.8603 -0.5146 1.0024 


-60° 0.4900 -0.8699 


30 “ 0.B717 


45 ° o.tu 4 


6o° O.ROPfi 




1 « 1 

. 

cut 5 a 

0.9963 

-90.0° 

0.9968 

-75.3° 

0.9984 

-6o,6° 

1.0005 

-45. 8° 

1.0024 

-30.9° 

1.0035 

-15.9° 

1.0036 

- 0.6° 

1.0029 

14.3° 

1.0013 

29.5° 

0.9996 

44,6° 

H.QQ7Q 

- « // 1 ^ 

■ 

RQ.R 0 
* “ 

0.9968" 

74.9° 

0.9963 

90.0° 



X 

y 

- 90 ° 

0,0000 

- 0.9932 

-750 

0.2540 

-0.9610 

- 6 o° 

0.4932 

-0.8659 

-45° 

0.70a 

- 0.7119 

-30° 

0.8659 

- 0.5070 

1 

J-* 

VJl 

0 

' 0.9708 

-0.2638 

0 ° 

1.0070 

-0.0000 

150 

0.9708 

0.2638 

30 ° 

O.S059 

0.5070 

45 ° 

0.7021 

0.7119 

6 o° 

' ' 

oJkhp 
~ - 

n.pjfca 

75 ° 

0.2540 

0,9610 

90 ° 

0.0000 

0.9932 









0.3183 0.3033 



0.0000 0.0000 


. I 

15.8° 0.1633 0.1612 


30.3'' 0.3183 0.3033 


45.4° 0.4557 0,4i46 


I g£UA i 


75,2° 0.6241 0,5294 


90,0° 0.6595 0.5505 






























































































































-90° 

O.OOOp 

-0.9923 

0.993 

=-90.0° 

0.5192 

0.4608 

-75° 

0.2467 

-0.9614 

0.9925 

-75.6° 

0.4972 

0.4452 

-6o° 

0.4856 

-0.8696 

0.9960 

-60.8° 

0.4331 

0.3974 

-1*5° 

0.6935 

-0.72® 

; 

0.999® 

-46.1° 

0.3326 

0.3156 

-30® 

0.8588 

-0,5196 

1.0Q37 

-31.2° 

0,2031 

0.1990 

-15° 

0.9672 

-0.2792 

1.0067 

-16.1° 

0.0539 

0.0538 

0° 

1.0079 

-0.0160 

1.0080 

-0.9° 

0.1094 

0.108J 

150 

0.9756 

0.2496 

1.0070 

14.3° 

0.2773 

0.2672 

7/V) 

S'* 

a nsn 

n lin^n 

WiT//t 

1 ftftln 

*»WVTA 

on C° 

ft kin 7 

V# ■ 

ft lisdin 

V* TV 'W 

-J 

U50 

0.7095 

0.7049 

1x001 

44.8° 

0.5913 

0.5090 

6o° 

0.4990 

0.5622 

0.9962 

59.9° 

0.7072 

153 

75° 

0.2572 

0.9594 

0.9933 

75.0° 

0.7926 

0.6212 

90° 

0.0000 

0.9923 

0.9923 

90.0° 

0.8204 

0.6342 



wns rm - telociw mbtoibtoiok ojt pbcstiz [m^- 0 . 3 ; a - ~eo°] 


A 

T 

m 

t « i 
1 * 1 

*■— ' 0 - 

f * 

1 

AA 

-90o 

0.0000 

-0.9922 

0.9922 

-90.0® 

0.3976 

0.3694 

-75° 

0.2452 

-O.9626 

0.9933 

-75-7® 

0.3771 

0.3529 

-&QO 

0.47*4 

-0,8741 

n acfUi 

70 

ft 7VT7 
"'■■/"'l j 

0.3024 

J150 

0.6859 

-0.7287 

1,0007 

-46.7® 

0.2218 

0.2165 

O 

O 

V 

0.8526 

-0.5317 

1.0048 

-31.9® 

O.0966 

0.0962 

-15° 

0.9637 

-0.2937 

1.0075 

-17.0® 

0.0510 

0.0509 

0® 

1.0078 

-0.0314 

1.0083 

-1.8® 

0.2179 

O.2129 

15® 

0.9791 

0.2349 

1.0069 

13.5® 

0.4o45 

0.3750 

30® 

0.8794 

0.4839 

1=0037 

PS, go 

0.5724 

0 s 4q68 

450 

0.7173 

0.6965 

0.9998 

44.2® 

0.7393 

0.5945 

< 5 o® 

0.5060 

■m 

0.9960 

59.5® 

0.8785 

O.6600 


o.a6i4 

0.9582 

0.9932 

74.7® 

0.9714 

O.6968 

90® 

0.0000 

0.9922 

0.9922 

90.O® 

1.0041 

0.7086 
































































































































































1 


i 


1 


g i.~ vmcnr diotribotioh oh raoraxa [m^*. 0.5; a * -eo°] 


1 z 1 arg b 


0.0000 

- 0.9795 

0.9795 

-90.0® 

0.7670 

0.2241 

- 0.9559 

0.9518 

-76.8® 

0.6575 

0.4368 

- 0 ,SB 4 o 

0.9561 

- . . i 

-63.7° 

0.5361 

0.6459 

- 0.7614 

0,9984 ' 

-49.7® 

0,3771 

0.8229 

-0.5256 

1.0100 

-35.4® 

0.1673 

0.9553 

- 0.3557 

1.0204 

-ao. 6 ® 

0.0919 

1.0219 

- 0.0920 

1.0260 

- 5 . 1 ® 

0.4089 

1.0069 

0.1903 

1.0247 

10 . 7 ° 

0 . 8025 - 

0.9077 

0.4576 

1.0165 

26.8® 

1.2953 

0.7372 

0.6819 

1.0042 

42 . 8 ® 

1.5824 

0.5162 

0.8466 

, 0.9917 

58 .e* 

2.5275 

0.2649 

0.9463 

0.9827 

M |. l._ 

3.0256 

A /WAA 
UlWW 

a r»*»nc 

a mnc 

voi^ 

flA AO 

1 

7 5C/YI 


0.9299 


0.9495 


tirt *ph Ha vow 1 
































































































































iabie n.- m'xxrriQX OF mcmix jbcm cutcnc m WCTimT ^ ^ r _, r- ^ 

[5^- 0,6; a * -20°J $#©3 XU.- YKLOCEKf IOSTOnaOTIOff 0® JKBTHI [l^s 0*7; a » 0°J 


A 

Z 

y 

H 1 

Mg S 


> 

X 

y 

1 1 1 

Mg S 

9 

" | 
AA 

-90o 

0.0000 

-0.9715 

0.9715 

-90.0® 

-90® 

0.0000 

-0.9494 

0.9494 

-90.0° 

2.7456 

0.9429 

-75® 

0.3395 

-0.9501 

0.9729 

-77.6® 


O 

IT \ 

I 1 — 

I 

0.2243 

-0.9281 

0.953** 

-76.4® 

2. 6081 

0,9337 

-6o® 

0.4-129 

-0.68SS 

0^9801 

-65.1® 

— 

-6o° 

0.4475 

-0.8618 

0.97U 

-62«6® 

2.2270 

0.9121 

*450 

0.6109 

-0.7865 

0.9959 

-52.2° 


JJ50 

0.6644 

-0.7420 

0.9960 

-48,2® 

1.6867 

0.S602 

O 

O 

V 

0.7997 

-0.6302 

1.0182 

-33.2° 

— 


-30® 

0.8606 

-0.5579 

1.0256 

-33.0® 

1.0916 

0.7374 

1 

•65 

0 

Q.9U70 

-0.4122 

1.0328 

— 

-23.5° 

-15® 

1,0063 

-0.3043 

1.0513 

-16.8® 

0.5127 

0.4562 

00 

1.0317 

-0.1432 

i.o4i6 

- 7.9° 

0® 

1.0620 

0.0000 

1.0620 

0.0® 

0.0000 

0,0000 

150 

I.O296 

0.1506 

i.o4o6 

8.3® 

15° 

1.0063 

0.3043 

1.0513 

16.8® 

0.5127 

0.4562 

30° 

0.9335 

0.433S 

1.0294 

24.9® 

30® 

0.S606 

0.5579 

1.0256 

33.0® 

1.0916 

0.7374 

1*50 

0.7565 

0.6701 

1.0106 

4i.5® 

45® 

0.6644 

0.7420 

0.9960 

4 *. 2 o 

1.6867 

0.8602 

6o« 

0.5263 

0.8392 

0.9906 

57.9® 


60® 

0.4475 

0.s6l8 

0.S7H 

62.6® 

2,2270 

0.9121 

75® 

0.2639 

0.9385 

0.9763 

74.0® 

75° 

0.2243 

0.9281 

0.9584 

76.4° 

2.6081 

HjS| 

90* 

0.0000 

0.9715 

0.9715 

90.0® 


90° 

0.0000 

0.9494 

9,9494 

90.0® 

2.7**56 



0.9429 


u> 

rs> 


f 


* 


< 
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> 

mnw 

mill 

HI 

arg * 

-90® 

0.0000 

-0.9499 

0.9499 

-90.0° 

-75° 

0.2151 

-0.9306 

0.9551 

-77*0° 

-60° 

0.4301 

-0.8693 

0.9699 

-63,7° 

-45° 

0.6409 

-0.7591 

0.9935 

-49.8® 

-30° 

0,8358 

-0.5882 

1*0220 

-35.1° 

-15* 

0.9897 

-0.3482 

l.o4gi 

-19.4® 

o« 

1.0607 

-0.0512 

1.0619 

- 2.8® 

1 

15° 

I.0227 

0.2582 

1,0586 

14,2° 

30° 

0.8860 

0.5262 

1,0300 

30.7° 

450 

0.6893 

0.7*3 

0.9998 

46 > 

6o° 

0.4663 

0.8545 

0.9732 

61.4® 

750 

0.2341 

0.9270 

0.9561 

75.8° 

90° 

0.0000 

0.9499 

0.9499 

90.0° 


} 


itf bib rev.- .VEioorrr DiamiBOTiOff os mams [m ot - 0.7; a - -io°] 


AA 


0.8236 


0.8184 




1.0373 


0.9115 


0.7154 


o.4s64 


0.2447 


0.0000 


1*1 

arg z 

0.9519 

-90.0° 

0.9565 

-77.5° 

1 ft ft7ftO 

Va 

-0 1 “»e 

0.9921 

-51> 

1.0195 

-37 .2* 

1.0463 

-21.9° 

1.0622 

- 5.5° 

1.0585 

11.5° 

1.0360 

28,4° 

.i.oo4g 

44.6° 

0.9769 

60.1° 

0.9584 

75.2° 

0.9519 

90.0® 




































































































































34 


HACA TS Vo, 1171 


TABLE XV.- DISTORTION OF PROFILE FROM CIRC EE BT TS TEN’S METHOD [m^- 0.7 ] 


A 

X 

y 

i * 1 

■BHI 

- 90 ® 

0,0000 

-1.1541 

1.1541 

-90.0® 

- 75 ® 

0.1445 

-1.1347 

1.1439 

- 82 . 7 ° 

- 6 o® 

0.3075 

-1.0661 

1.1097 

- 73 . 9 ° 

>450 

0.4892 

-0.9250 

1.0464 

-62.1® 

e 

O 
t«"\ 
. t 

0.6659 

-0.6926 

0.9608 

- 46 . 1 ® 

- 15 ° 

0.7971 

-0.3731 

0.8801 

- 25 . 1 ° 

o« 

0.8459 

0.0000 

0.8459 

0 . 0 ® 

150 

0.7971 

0.3731 

0.8801 

25 . 1 ° 

0 

0 

0.6659 

0.6926 

0.9608 

46 . 1 ® 

450 

0.4892 

0.9250 

1.0464 

62.1® 

60® 

0.3075 

1.0661 

1.1097 

73 . 9 ° 

75 ° 

0.1445 

1.1347 

1.1439 

82.7® 

90® 

6.0000 

i. 154 i 

1.1541 

90.0® 
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rig*. i,s 



Tigure 2.- Shape of profile. U “ 0.1; a - -20° 



HMM TH HO. 1171 Tige. 3.4 
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fig*. 5,6 




Figure 5.- Sh«pe of profile. » 0.7, 
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fig*. 7,8 
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ARC X 


Figure 7.- Maoh number as a function of the argument s. Hjkj » 0.7. 



Figure 8 U&oh number as a function of the argument x. il m » 0.3. 
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rigs. 9,10 




figure 10.- Shape of profile. 1 £qo « 0.7; a- 0°. 





